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We explore the advantages offered by twin light beams produced in parametric down-conversion for precision 
measurement. The symmetry of these bipartite quantum states, even under losses, suggests that monitoring 
correlations between the divergent beams permits a high-precision inference of any symmetry-breaking effect, 
e.g. fiber birefringence. We show that the quantity of entanglement is not the key feature for such an instrument. 
In a lossless setting, scaling of precision at the ultimate 'Heisenberg' limit is possible with photon counting 
alone. Even as photon losses approach 100% the precision is shot-noise limited, and we identify the crossover 
point between quantum and classical precision as a function of detected flux. The predicted hypersensitivity is 
demonstrated with a Bayesian simulation. 

PACS numbers: 42.50.-p.42.50.St,06.20.Dk 



I. INTRODUCTION 

We wish to promote the fitness of two spin-j systems com- 
bined in an overall spin zero singlet as a resource for metrol- 
ogy. A photonic implementation of such states is readily 
and scalably generated using stimulated parametric down- 
conversion (PDC) [ 1 1. Both their rotational symmetry J21 and 
persistence of entanglement ||3] |4] under photonic loss chan- 
nels has made the singlets natural candidates for quantum key 
distribution [5 1. In this letter we further highlight their utility, 
this time in a parameter-estimation protocol. Generally, such 
protocols are rated by the precision (or uncertainty) associated 
with unbiased estimation of the unknown parameter, and how 
quickly this precision is lost under relevant decoherence. We 
will show that given ideal conditions the singlets allow a pre- 
cision scaling at the Heisenberg limit (the ultimate limit for 
linear quantum processes, and for which noise scales as 1/iV 
with respect to the light intensity or particle number N). Un- 
der incoherent photon loss measurement precision is naturally 
degraded, but at a much gentler rate than other proposals [6| 
where the decay can be exponential in N. (Recently, the role 
of photon losses in optical precision experiments was exam- 
ined carefully Q.) 

Consider an instrument broken down into three components 
[8|: pure probe state \j), unitary evolution U = exp(— i<pHt) 
under a time-independent hermitian Hamiltonian H, and com- 
plete projective measurements, M = 2i m iNX*l> where 
= §ij. We wish to infer evolution time t from frequen- 
cies of outcomes rrij, but t may equally represent an interfer- 
ometer phase, a magnetic or gravitational field, or some other 
real-valued continuous variable. Extrapolating from a set of 
parameter-dependent measurements to an estimation of that 
parameter can be a challenging task. The probability distri- 
butions for individual outcomes are often non-Gaussian, hav- 
ing multiple peaks or broad tails. In addition to designing 
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a good estimator from measurement data, it is as important 
to employ measurements sensitive to small changes in the 
unknown parameter. An optimal measurement should also, 
if possible, achieve highest estimation precision for all pa- 
rameter values. Conditions for identifying optimal measure- 
ments were identified in early work on quantum Fisher infor- 
mation. One parameter-independent approach uses canonical 
measurements, but these are hard to implement directly J5). 
In other proposals, measurements are optimal near particular 
'sweet spots' in parameter space, and require multi-step adap- 
tive measurements to exploit them iflOll . 

We first present a protocol in a decoherence-free setting 
for which our proposed measurements are optimal, parameter- 
free and practicable; they may be realized in a laboratory by 
photon-counting or spin projections. We extend the analysis 
to consider realistic decoherence. Translating our protocols 
to a quantum optics setting, characterized by a PDC source, 
we evaluate the effect of photon losses (in transmission and 
detection) on precision. It emerges that spin-projection mea- 
surements are no longer optimal or parameter independent. 
Nonetheless, precision is always at least as good as the classi- 
cal upper bound (noise scaling oc iV~ 1//2 , called the shot-noise 
limit) for any loss. 



II. INPUT STATES 

Let us introduce the state in the spin representation ifTTl . 
Our probe is a bipartite maximally-entangled spin singlet: 

|^'>} = — L_ J2 (-l) j - m \j,m) za ®\j-m) zb (1) 

m=—j 

labeling component spaces 'a' and '6'. The state has total 
spin J = 0: {ip ( 3) \J 2 \^ 3) ) = where J = ( J a + J b ). It 

is 'rotationally' invariant |2| under Ua^ (<?i) <8> (32) when 
.9i = 92 € 5(7(2); its description in Eq. ([T| is identical in any 
spin basis, e.g. z n- x M> y, and it has properties that change 
only with the relative transformation l a ® (9i)U^ {92)- 
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FIG. 1: Instrument composed of PDC source [ Q, linear-optical el- 
ements (mirrors, phase element </>, 50/50 beam-splitters l bs', po- 
larising beam-splitters 'pbs\ half-wave plate 'A/2'), and photo- 
counters. In microscopy or interferometry, modes a may be con- 
sidered an ancilla/reference, only the photons in modes b inter- 
act with the 'sample'. The down-conversion Hamiltonian is H = 
^(a^bl — albf^ + h.c. with k a coupling strength for the nonlinear 
process, producing light in four optic modes via application to the 
vacuum: exp(— iHt)\0). Given effective interaction time r = «i, 
and disregarding losses, the source produces a mixture of spin-j sin- 
glets, Eq.^TJ, with weighting (2j + 1) tanh 4j r/ cosh 4 r, and each 
singlet has a photon number N = 4j. Overall, the intensity is 
(N) = 4sinh 2 r and A 2 iV = cosh4r — 1. Weightings of higher- 
photon-number singlets increase with r but the large A 2 N indi- 
cates a severe flattening of the distribution. We note the immunity 
of counting measurements to undesired path differences occurring 
outside of the 'interferometer' block. For fiber calibration the ap- 
paratus would be simpler, comprising polarizing beam-splitters and 
photon detectors terminating both test and reference fibers coupled 
to the PDC source (without an explicit interferometer or A/ 2 plate). 



There is an application here for relative measurements made 
between non-local observers; global phenomena are excluded. 



IV. QUANTIFYING PRECISION 



To establish this scheme's performance we employ classical 
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Fisher information I c i ((f)) = J2a b Pa,b 
provides a distance metric in probability space, e.g for e<C 1 
the distance between Pab(4 i ) an d Pab(4> + e ) is e VPi- A 
lower bound on the variance of any unbiased estimator of <f> is 
Id 1 G2- Gi yen an unknown rotation by <f> about the 'by' axis 
a short calculation (Appendix A) using Eq.|2]) gives I c i (<f) = 
4j (j + 1) /3, independent of tf> and achieving the quadratic 
scaling characteristic of the Heisenberg limit (identifying a 
spin-j state as a composite of 2j spin one-half particles, the 
singlet then has particle number N — 4j and I c [ oc N 2 ). 

The quantum Fisher information I qu , provides a saturable 
upper bound, I qu > I c i. It is a function only of probe and 
dynamics, assuming the best possible measurement without 
defining it explicitly. For pure probe states, I qu /A = A 2 H = 
(H 2 )-(H) 2 ED. Using again H ^ J yb , I qu = 4(J 2 



4(J?>/3 



ybl 



4j(j + l)/3 by symmetry. Remarkably I qu = I C . 



and our original measurement choice is optimal and indepen- 
dent of the value of the unknown parameter 4>, a preferred 
property of any parameter estimation scheme JSJ [14). We 
make at this point some observations about entanglement. 
First, making a comparison between the probe | V'o"' ^ ) an d Bell 
states IV^ 2 ') numbering 2j (so total particle number 4j is 
the same), we see that the former has much less entangle- 
ment than the latter, log 2 (2j + 1) e-bits versus 2j e-bits, but a 
greater value for I qu , 4j(j+l)/3 versus 2j. Second, although 
the Hamiltonian only operates on the b modes, no phase in- 
formation is imprinted locally as the reduced state /?w is al- 
ways maximally mixed - nothing is learned by measurements 
exclusive to the space in which dynamics occur. One might 
believe that the dependence of precision on bipartite measure- 
ments is due to entanglement. However, we will show it is 
retained under a disentangling decoherence channel - exhibit- 
ing non-locality without entanglement. 



III. PHASE-ESTIMATION PROTOCOL 



V. SIMULATION 



Examine now the estimation of a relative phase 4> ac- 
cumulating between a and b. Breaking the symmetry of 
\tp0^) by applying 1 <g> uf?\ yields the state \i>^ ) ) = 

cxp (—i(/)jyb) IV'o^) if we choose rotation about y. In the 



paradigm above, H H> jyj, and ( h> ^. Then make projec- 
tive measurements onto the basis of z-eigenstates for a and 
b, i.e. M i— !► j za ® J z i,. The probabilities for the (2j + l) 2 
measurement outcomes are 



P AB (<)>) = (j,AL a (j,B\ b \^) 



2 d 



B,-A 



(2j + 1) 



(2) 



Note that Pab (</>) — Pab (~<P)> an d to eliminate ambiguity 
the range of <f> can be restricted to [0, ir]. 



The quantum bound (8(j)) 2 > I q ^ on mean-squared er- 
ror is attainable given infinite repetitions of the experiment 
(using maximum-likelihood estimation), but superior conver- 
gence of the high-j singlets may also be shown when the 
number of samplings is relatively small, the pre-asymptotic 
regime. A Bayesian protocol, considered in Fig. [2] com- 
pares performance of higher- j singlets with Bell singlets, so 
that the same total energy is detected in both cases. We 
use data sequences with k updates of the (uniform) prior 
distribution Pq = 1/n given measurements {Ai,E>i} — 
{Ai, B 1: A 2 , P> 2 , ■ ■ ■ , Ak, -Bfc}. Each sequence leads to a sin- 
gle <f> estimate that is the mean of the final posterior distribu- 
tion Pfc. The conditional probability P({Ai, Bi}\<f>) is deter- 
mined from Eq. 0. The update rule is Pu (4>\{Ai, B{\) oc 
P({A i7 Bi}\(f))PQ{4>) and measurement events are indepen- 
dent, so P(A l ,B 1 ,A 2 ,B 2 \ct>) = P(Ai,B l \<f>)P(A 2 ,B 2 \4>). 
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FIG. 2: Simulation with performance of spin-j singlets (red/above) 
juxtaposed with spin-1/2 Bell pairs (blue/below). Like coherent state 
light, the Bell states have a precision I qu that scales linearly with 
particle number. There were 250 independent phase estimates con- 
tributing data points to each of the six charts. A single estimate is 
obtained after a sequence of k = 20 Bayesian updates for probes 
with j = 2,3,4 using the distribution of Eq.{2|. The true value of 
<f> is indicated by black dashed lines. Maintaining the same par- 
ticle number resource for the Bell probes as the higher-j singlets 
equates with keeping jk constant. Experimental error bars approach 
the theoretical limit 8<f) m — 1/ y/kl qu as particle number lOOOjfc 
increases; red bars are smaller than corresponding blue ones (below 
them) indicating superior precision of the higher-spin singlets in this 
pre-asymptotic regime. 



VI. LINEAR OPTICS AND PARAMETRIC 
DOWN-CONVERSION 



Now we translate our previous arguments into an optical 
context using an isomorphism between spins and a pair of har- 
monic oscillators (Schwinger representation). As such, one 
application that suggests itself is in measuring birefringence 
of optic fibers to high precision by comparing a test sample 
with a calibrated reference. The PDC process produces entan- 
gled photons by interaction of a pump laser field with a non- 
linear birefringent crystal. The output is shared among four 
optic modes, labeled {a^, a v ,bh,b v } for spatial directions a, b 
and horizontal/vertical polarizations h, v. (We use the same 
notation to denote the associated bosonic annihilation opera- 
tors.) Applying the Schwinger representation to spatial mode 
a: J +a = a[a v , J_ a = a h a\, and 2J Z = (a} h a h -a\a v ) 
(the number difference). Also J 2 = (n a /2)(n a /2 + l) where 
h a = (a\ L ah+a\,a v ). Spin quantum numbers map onto pho- 
ton numbers as 2j a = {n ha + n va ) and 2m a = {n ha -n va ). 
We can now identify each of the elements of our idealized 
parameter-estimation protocol with a realistic optical source, 
as explained in Fig. [T] Measurement data can be grouped 
by photon counting according to values of j a ^, equivalent 

to post-selection onto the space of a particular \i[>q^). Op- 
timal correlation measurements J za g3 J z t, are reconstructed 
in each j a b-labeled subspace, also from photon counting be- 
cause 2 J za = hha — n va . We acknowledge here that efficient 
photon counting is generally a non-trivial task lfl~5ll . 



VII. LOSSES 

A realistic analysis must incorporate relevant decoherence; 
for optic processes incorporating photon counting the impor- 
tant mechanism is that of photon loss, in transmission and de- 
tection |4|. Both loss types are effectively modeled by plac- 
ing partial transmission (r/ < 1) beam-splitters in the four 
optic modes in front of perfect detectors, splitting incoming 
photons into two output modes: the mode transformation is 
a i y ^Jr)a + y/T— r/e where e is an annihilation operator for 
the ancillary loss mode. The photons syphoned out of the 
transmission mode in this way are then traced over. Referring 
either to modes a or to modes b, if losses are polarization- 
insensitive (i]h=Vv), the l° ss channel C^<E)C^ commutes with 
any U E SU{2) on the same spatial path: t\ ® O^UpW] = 
U(C^ <X> Cl\p])W for any p. This has two important con- 
sequences. Firstly, the non-unitary decoherence due to loss 
and the unitary <f> rotation in mode b may be treated indepen- 
dently with impunity. Secondly, after losses each component 
of the PDC state, in a subspace labeled by (j a ,jb)> retains its 
symmetry under transformations Ua'(g) ® U^ b \g). This 
implies a simple, block-diagonal structure for the mixed lossy 
state: 

p V*M = l Jj.ji,) sl (j*j>) 4 (3) 

J=\ja-h\ 

The n < j a ' j ^ e [0, 1] are weighting factors H and O^'^ are 
density operators proportional to identities in each (2 J + 1)- 
dimensional orthogonal subspace labeled by total spin J. As 
symmetry is preserved under loss, the state of Eq. Q retains a 
suitability for relative measurements between spatially sepa- 
rated observers a and b. We stress that for imperfect transmis- 
sion rj < 1 there is non-zero occupancy probability for spaces 
labelled (j a ,jb) where j a ^ jb- See Fig.|3ji). Measurements 
with n a 7^ rib need not be discarded, they also contribute to 
overall precision. 

VIII. COMBINING SUBSPACES AND ULTIMATE 
PRECISION 

By post-selecting onto specific photon numbers (after 
losses) and filtering the data sets we can focus on a particular 
(j a ,jb) subspace and analyze its contribution to overall preci- 
sion. Fortunately, both I qu and I c i are additive, so total Fisher 
Information per measurement is the average of the contribu- 
tions in each subspace weighted by the probabilities of post- 
selecting each subspace. For the lossless case with N = n a +hi, 
we know I qu = 4j(j + l)/3 (in a space with n a> b = 2j pho- 
tons) and the ensemble result is (I qu ) = (N 2 )/12+ (AO/3. 
There are associated gains in precision as (N) = 4?? sinh 2 r 
increases with r in the low loss regime. For more significant 
loss there is a precision trade-off as larger r are also associated 
with greater values of decoherence 9= (1 — rj) tanhr within 
each subspace. This may be viewed as higher-photon-number 
spaces being seeded initially, then as photons are lost these 
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populations make an incoherent contribution to the those of 
lower photon spaces, 'feeding' them from above with mixed 
state components (bad for precision). We illustrate the effect 
of decreasing transmission 77 on both subspace weightings in 
Fig.|3ji), and precision within a subspace, in Fig. [3jii) (explic- 
itly the j a —ji, — 2 subspace). 

Given a general mixed state p = J2 P Pp \p) (p\ Fisher Infor- 
mation for unitary evolution under H i-» J y h is 



(I VU )=2J2 



{Pq-Prf 
Pq + Pr 



{<l\Jyb\ 



(4) 



By observing that the PDC is Gaussian, and that loss channel 
and interferometer components act as Gaussian operations, 
this functional can be evaluated directly using phase-space 
methods |[T6ll . Specifically, we employ the fact squeezed light 
subject to incoherent photon loss is formally equivalent to a 
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FIG. 3: (i) Under photon loss the PDC state becomes 
Y,P(n a ,n b )p ( - ja ' jb \ where j a = n a /2,j b = n b /2, pWo.ii.) is given 
in Eq. {3}, and P(n a ,nt) is a weighting factor for the correspond- 
ing subspace |4|. The three plots show P(n a ,n b ) with interaction 
r ~ 1.83 and n a ,b < 16. As transmission 77 is reduced, the weight- 
ings of j a 7^ jb spaces increase and n a ,b 3> 1 become less proba- 
ble. In (ii) loss in transmission and detection affects the estimation 
of (j> in the (n a ,rib) i-> (4,4) subspace. The contribution to (I c i) 
is 7^' 4 'P(4,4). For no losses (77 = 1), p 1 - 2 ' 2 ' is the singlet state 
|^Q 2 ')('i/)Q 2 ^| and I^f'^ = 8, contributing to the quadratic precision 
scaling derived for the lossless case. As losses increase, the preci- 
sion degrades and I^f ,4 ' quickly becomes a function of both <j> and 
the decoherence parameter 8 — (1 — 77) tanhr. Blue shading indi- 
cates precision below the upper limit for four uncorrelated photons 
in a classical lossless two-mode interferometer, i.e. I c i < 4. For 
<f> « 0.5 supra-classical precision (red shading) is possible for deco- 
herence 9 < 0.2. This corresponds to loss < 26% for r = 1, an 
interaction value achieved in real experiments (3) . The maximum 
/$ 4 ' 4 ' for each 6 value with photon counting measurements is pro- 
jected on to the back wall as a grey silhouette - only slightly inferior 
to /gu' 4 ' , the best possible for any detection scheme (red dashed line). 



squeezing operator (with a modified parameter) applied to a 
thermal state. This allows for the a full identification of the 
spectrum for the lossy PDC state, and for Eq. (|4]) to be eval- 
uated directly by evaluating nonzero contributions to the sum 
(Appendix B). We arrive at an expression for the complete 
Fisher information in terms of transmission 77 and detection 
flux (N) = 4?/sinri 2 r: 



(N^)+4(N) {N)(4r,+(N) 

> \J-qu) — 



12 



> 



8+4(l-77)(AT) 8+4(N) 



(5) 



where the upper and lower limits correspond to 77 = 1,0 re- 
spectively. Therefore, even in the worst case, losses approach- 
ing 100%, any non-zero flux gives (I qu ) ~ (N) /4 and shot 
noise limited precision scaling 8<j) ~ (TV) -1 / 2 (in stark con- 
trast to the exponential deterioration in performance of other 
schemes [6|). As losses increase, flux (N) can be maintained 
by turning up the interaction r, but the Fisher information will 
inevitably deteriorate because of its separate dependence on r\ 
and (N) above. It is not known whether optimal measure- 
ments in the lossy case can be independent of the true pa- 
rameter value tf>, certainly photon correlation measurements 
Jza 8) J z b are no longer optimal, see Fig.^n). From Eq. |5]) 
it is seen that our scheme performs better than a sample illumi- 
nated with coherent light of the same detected flux, for which 
(I qu ) = (N)/2, when 77 > 1/2 + 1/(2 + (N)). Transmis- 
sion must be better than 50% in the high flux limit to obtain 
a quantum advantage with PDC. Fig. [4] shows the scaling of 
(Iq U ) with detected flux (N) when fitted locally to a power 
law: (I gu ) oc (N)\ 
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FIG. 4: The scaling of (I qu ) with the detected flux (N) is examined 
in the plot, for various values of the transmission parameter 77. A 
logarithmic scale is used for (N) on the horizontal axis. Scaling 
parameter 7 is defined by 7 = d \og({I qu )) / d\og({N)) , and locally 
(Iqu) °c (AO 7 - We see that for our protocol 7 varies between 1 (shot- 
noise scaling) and 2 (Heisenberg scaling), i.e. our scheme always 
scales better than the best classical interferometer using uncorrelated 
particles. For moderate losses 77 w 0.8 the optimal flux is between 1 
and 10 photons, well within the capabilities of some current number- 
resolving photodetection technologies 1 15 1. 
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IX. SUMMARY AND OUTLOOK 

We have presented a parameter-estimation protocol with 
several strengths. In a non-dissipative environment, 
Heisenberg-scaling is achieved with simple fixed measure- 
ments. Our scheme can be implemented using PDC and lin- 
ear optics, and under severe dissipation (approaching 100%) it 
is capable of precision scaling that approaches the shot-noise 
limit asymptotic ally from above. There exist a variety applica- 
tions, such as fiber calibration, the phase microscopy of frag- 
ile biological specimens, and optical gyroscopes for GPS-free 
navigation. We will, in a forthcoming work, extend this ap- 
proach beyond a single phase estimation to all of SU(2), pa- 
rameterized by three Euler angles, and utilizing multiple spin 
correlations. The increasing asymmetry of the state with the 
magnitude of one-sided rotations has a quantifiable utility in 
reference frame alignment 1 17 1. 
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Appendix A: Derivation of J c ; for pure-state probes and no 
dissipation 

Referring to the parameter-estimation paradigm described 
in the main text, we set H M> J y i, and t i-» <f>, and 
M H> J za (8 J z b. Under the unitary dynamics, the j-singlet 



probe state evolves to \4' l 
are (2j - 



0')\ 



exp 



There 



1) outcomes, and the probability distribution for 

2 

these is of the form P AB 0) = (j, A\ za (j, B\ z b |Vv-) 

d%_ A (0) 2 /(2j + l). 

The classical Fisher information is defined to be, 
!d (4>) = Y.A.B P A,B[d\n(P A . B ) /d(j)] 2 . The deriva- 



tives dPAB/d(j) are given by dg _ A (0) / (2j + 1) 



{-A) d»>_ A _ x (0) -N%> {-A) d^_ A+1 (0) 

ming over all detection outcomes gives I c i (</>) 
4j + /3, and this is independent of the value of < 
equal to I qu . 
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Sum- 
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Appendix B: (I qu ) for the PDC state in the case of equal loss 
rates for every mode 

Referring to the parameter-estimation paradigm described 
in the main text, we set H n- J b and t (f>. Now an op- 
timal measurement is assumed, but is not made explicit, and 



we derive the best precision possible in principle. Our probe 
state is assumed to be the PDC state corresponding to effec- 
tive interaction time r. Photon loss arising from transmission 
and detection is assumed to act independently on all four optic 
modes, and the same transmissivity parameter 77 is assumed in 
each case. 

The calculation of (I qu ) is presented in four stages: (i) the 
lossless PDC state is written in terms of one-mode squeezing 
operations; (ii) the action of the loss super-operator on the 
squeezed vacuum is explained; (iii) eigenvalues and eigenkets 
for the lossy PDC state are identified; (iv) (I qu ) is evaluated 
in terms of t, rj. 

( i) The lossless PDC state represented as a Gaussian state: 
The PDC state is of the form, 



00 

|^ pdc )(x^ % A^TTA"e^|C /2) ) ) 

where the phase factor cp is the phase of the pump field, 
A = tanh (t) corresponds to the strength of the non-linear 
interaction, and the normalized photonic n/2-singlet state is 

|V>0 n/2) ) OC Em=0 I" ~ m > m ' m ' U - m )a h a^b h K- 

Next, we rewrite the PDC state in terms of Gaussian op- 
erations acting on the vacuum, and specifically in terms of 
squeezing processes acting in single modes: 

l^pdc) = u sio Mu h a l b Ualb h 

oS ah (r)S av (-r)S bh (r) |0>&„ (-r) |0), 

where C/^ = cxp (7r/4) (jtfb — db^ denotes a 50:50 beam- 
splitter operation, J7 g i b a i(<y3) = exp(iNip/2) accounts for the 
global phase factor, and the squeezing operations are de- 



fined by S ah (t) = exp 



<r/2) (4 2 



S a „ (-r) 



exp [(t/2) (aj 2 — a 2 ,)] , and similarly for modes b. 

( if) The loss channel as a Gaussian operation and its action 

on the squeezed vacuum state: The loss super-operator for 

00 

mode c G {a h , a v ,b h ,b v } acts as C^[p] = L nP L n > with 

71=0 

L c n = (1 — rj) n / 2 c n if c / 2 /\fn\, where the loss in mode c is 
(1 — rj) G [0,1]. We note that £Jj is a Gaussian operation, 
and therefore that the action of £Jj on the squeezed vacuum 
state for one mode can readily be computed using standard 
continuous-variable methods. We make the following formal 
identification: 

£*[£« (±r) |0)] = S (±r eff ) p th (N)& (±r cff ) , 

where /5th denotes the thermal mixed state p t h oc 
Yl^Lo X n \ n )( n \> m tne Fock basis, and x = N/ (l + N). 
The thermal intensity parameter N and the effective squeez- 
ing parameter r e g are related to the PDC and transmission 
parameters as follows: 



Teff = 



N = 



= im 



4 
-1 



p 



2t 



1 — 77, and M = rje' 



P 
M 
y/PM 

2r + 1 



6 



( Hi) Eigenvalues and eigenvectors of the lossy PDC state: 
Since the photon losses that arise due to transmission and de- 
tection are assumed to be polarization insensitive, they com- 
mute with the interferometric stages of our protocol, and can 
therefore be assumed to act exclusively at the PDC source. 
Furthermore, because the same value for the transmission pa- 
rameter 77 is assumed for all four modes, the PDC state can be 
written as follows, after accounting for all losses: 



; lossy -f7cff [p th (Nf 



PpBC : 



U, 



eff' 



where, 



U cS = 



°S ah (r cff ) S av (— T cff ) Sb h (r e ff) Sb„ 



(-Teff) 



Definitions for U hs , £/ g i bai and S are as in subsection (i). 
Parameters T c ff and N are as in subsection (ii). 

The eigenkets of Pp^c are therefore given by 
Ucs\Ah,A v , Bh,B v ), for every product of Fock states 
in four modes: \A h ,A v ,B h , B v ) ahavbhbv . The corresponding 
eigenvalues are of the form x A h+A v +B h +B v j + . 

(iv) Evaluating (I qu ) for the lossy PDC state: For a 
mixed probe state p ((f)), subject to a unitary evolution 

exp (— i(j}H\ , the quantum Fisher information is given by 

(I qu ) = tr [(dp/ defy [dp /deft}, where RJ 1 [dp/d^] de- 
notes the symmetric logarithmic derivative of p. The value 
of (Iqu) is independent of <f>, and we set <p — 0. By diago- 
nalizing the probe state, p = J2 P Py \p) (p|> the quantum fisher 
information can be evaluated using the expression, 



q,r 



(Pq - Prf 
Pq + Pr 



[q\H\r) 



In what follows H (£>u&1 — "V>h) /2*, corresponding to 



the Hamiltonian for a y-axis rotation in modes b, while the 
eigenkets and eigenvalues are those, given in subsection (iii), 
for the PDC state after the loss super-operator has acted on all 
four modes. 

To evaluate (I qu ), we treat it in stages. There are eight sum- 
ming indices, Ah,A v ,Bh,B v and A' h ,A' v ,B' h ,B' v , correspond- 



ing to possible photon numbers in the four modes. Next, 



(pg - Pr) 2 

Pq + Pr 



(i + nY 



A h +A v +B h +B v _ x A' h +A' v +B' h + B' v 



I + i„+B h +B v + x A' h +A' v +B' h +B' v 



By acting U c s on the Hamiltonian, the final component of the 
sum is found to be given by, 



(q\H\r) 



1 ? 



1 

16 



sinh (2tpi 



where T — 



x (A' h A' v B' h B' v \T\A h A v B h B v ) 

hd v + b v b h + a\dX — b\b\, and terms 
which do not contribute to the sum have been dropped. 
For a given value of the summing indices, at most 
one component of T can generate a nonzero con- 



tribution. Hence, 



(A' h A' v B' h B' v \f\A h A v B h B v ) 



three similar terms. 



\(A' h A' v B' h B' v \a h a v \A h A v B h B v )\ 
Furthermore, each of these four terms provide an equal total 
contribution to the sum. 

Now (Iqu) can be evaluated directly: 

_ sinh 2 (2r eff ) (2N + 1) 2 
{ qu) ~ 2 (1 + 2N + 2iV 2 ) 

[77 cosh 2 (t)] [2rj sinh 2 (r)] 
1 + (1 — 77) [2?ysinh 2 (r)] 

Finally, it can be recast in terms of the total detected flux 

(N) = 477 sinh 2 r, 



(Iqu) 



(N) (4 V +(N) 
8+4(l- v )(N) 



Using this result, the scaling of (I qu ) as a function of the 
detected flux and the transmission parameter is examined in 
Figg] 
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